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g : ABSTRACT 
O 

From the peak of a gravitational microlensing high-magnification event in the 
A component of QSO 2237+0305, which was accurately monitored by the GLITP 
collaboration, we derived new information on the nature and size of the optical 
l^-band and i?-band sources in the far quasar. If the microlensing peak is caused 
by a microcaustic crossing, we firstly obtained that the standard accretion disk 
^ ' is a scenario more reliable/feasible than other usual axially symmetric models. 

^ ■ Moreover, the standard scenario fits both the V^-band and i?-band observations 

■ with reduced chi-square values very close to one. Taking into account all these 
^ ■ results, a standard accretion disk around a supermassive black hole is a good can- 
O ■ didate to be the optical continuum main source in QSO 2237+0305. Secondly, 

• using the standard source model and a robust upper limit on the transverse galac- 

tic velocity, we inferred that 90 per cent of the l^-band and i?-band luminosities 
2 ■ are emitted from a region with radial size less than 1.2 10~^ pc (= 3.7 10^^ cm, 

c/3 . at 2 cr confidence level). 

> ' 
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1. Introduction 

In the optical continuum, QSO 2237+0305 is a gravitational mirage that consists of 
four compact components (A-D) round the nucleus of the deflector (lens galaxy). The light 
bundles corresponding to the components are passing through the bulge of the lens galaxy, 
and thus, if the galactic mass at the QSO image positions is mainly in stars, the optical depths 
to microlensing are as high as ~ 0.5 (e.g., Schmidt, Webster & Lewis 1998). In this scenario 
with large normalized surface mass densities, given a component, microlensing violent events 
will result from the source cither crossing a microcaustic, passing close to a microcusp, or 
traveling through a network of microcaustics (e.g., Schneider, Ehlers & Falco 1992). On the 
other hand, a microlensing violent episode in the light curve of a component can be easily 
distinguished from a intrinsic variation, since the intrinsic variability must be observed in 
all four components of the system with extremely short time delays (e.g., Wambsganss & 
Paczynski 1994; Chae, Turnshek & Khersonsky 1998; Schmidt, Webster & Lewis 1998). 
However, when the four brightness records of the QSO images have different non-flat shapes, 
a direct separation between the true microlensing fluctuations and the possible intrinsic 
variation cannot be achieved. In the case of four incoherent and non-flat observational 
trends, to find the true microlensing behaviours we must do some hypothesis on the intrinsic 
variability. 

Irwin et al. (1989) discovered microlensing variability in the quadruple system QSO 
2237+0305, and that first evidence was confirmed by other observers (Corrigan et al. 1991; 
Ostensen et al. 1996; Wozniak ct al. 2000a,b; Schmidt et al. 2001; Alcalde et al. 2002). 
In very recent years, two gravitational microlensing high-magnification events (HMEs) were 
clearly detected by the OGLE collaboration (Wozniak et al. 2000b) and corroborated by 
the GLITP monitoring (Alcalde et al. 2002). As each individual HME is directly related to 
the intrinsic surface brightness of the source, the two HMEs in the V band reported by the 
OGLE team were used to obtain two measurements of both the optical continuum (T^-band) 
source size and the rate of brightness dechne (Shalyapin 2001). Shalyapin (2001) compared 
the observational data included in each HME with the time evolution expected from an 
axially symmetric source crossing a single straight fold caustic. To describe the brightness 
distribution of the y-band source, he used a power-law model: /y(r) = 2^^/y(l+r^/i?y)~^^, 
which is determined by a typical radius Ry, a typical intensity ly = Iv{Rv)- and a power- 
law index py. The expected microlensing light curves depend on five free parameters, and 
the most relevant ones are At = Rv/V±_ and pv- We note that a direct measurement of the 
typical radius Ry is not possible, however, using some upper limit on the quasar velocity 
perpendicular to the caustic line (Vj_), it can be obtained a very interesting constraint on 
the source size as measured by means of the typical radius of the 2D brightness distribution. 
A crossing time of At 90 days is inferred from the HME observed in the light curve of the 



-3- 



component A, whereas a shorter crossing time of about 30 days is consistent with the HME 
corresponding to the image C. However, the whole HME of image A (from day 1200 to day 
1800) seems to be caused by a complex magnification, which is different to the single straight 
fold caustic magnification law. In other words, when the main portion of the source is far 
away from the fold caustic of interest (before day 1400 and after day 1600), there is evidence 
for a rare behaviour, and so, the fit to the whole microlensing event in the brightness record 
of the component A could give biased estimates of the parameters. Taking into account 
this perspective, we only take the results based on the HME corresponding to the image C 
as non-biased parameter estimates. The value of At 30 days together with the velocity 
constraint by Wyithe, Webster & Turner (1999) give an upper hmit on the F-band typical 
radius of Ry < 3 10~^ pc. A similar conclusion was obtained by Yonchara (2001), who 
analyzed the same microlensing event but using a different picture and a typical microlens 
mass of ^ 0.1 Mq (Wyithe, Webster & Turner 2000). On the other hand, Shalyapin (2001) 
found that the value of the power-law index is close to the validity limit of the model {pv ^ 
1) and the best-fit reduced is significantly greater than 1. 

Apart from the very recent papers by Shalyapin (2001) and Yonehara (2001), other 
previous works also discussed the size of the optical continuum source (e.g., Wambsganss, 
Paczynski & Schneider 1990; Webster et al. 1991; Wyithe ct al. 2000). These first studies 
are based on a poorly sampled HME which was observed in Q2237-I-0305A during late 1988. 
While the continuum mostly arises from a compact source, the line emission comes from a 
much larger region. Prom two-dimensional spectroscopy of QSO 2237-1-0305, Mediavilla et 
al. (1998) found an arc of extended C III]A1909 emission that connects the components A, 
D, and B. The observed arc is consistent with a source radius larger than 100 pc. 

The GLITP (Gravitational Lenses International Time Project) collaboration has moni- 
tored QSO 2237+0305 in the period ranging from 1999 October to 2000 Pebruary (Alcalde et 
al. 2002). The GLITP/PSPphotll photometry in the V and R bands (see Pig. 1 in Alcalde 
et al. 2002), showed a peak in the fiux of the component A and a relatively important gra- 
dient in the flux of the component C. These features are related to the two HMEs that were 
discovered by the OGLE team. The GLITP light curve for image A traced the peak of the 
corresponding HME, i.e., the maximum and its surroundings, with an unprecedented quality. 
Por example, the OGLE collaboration sampled the F-band peak at 19 dates, whereas the 
GLITP record included measurements of the F-flux at 52 dates. Moreover, the rest of global 
behaviours (components B-D) were accurately drawn from the GLITP photometry (in this 
paper, we will use the PSPphotll variant). The VR light curves of the four components A-D 
have been also analyzed from a phenomenological point of view, and the global flat shape 
for the light curve of Q2237+0305D suggested that the microlensing signal in D and the 
intrinsic signal are both globally stationary. In consequence of this result, it seems that the 



-4- 



global variabilities in A-C are unambiguously caused by microlensing. 

As mentioned here above, the whole HME of Q2237+0305A seems to be originated by 
a complex magnification law. However, in principle, the peak of the HME (just when the 
main portion of the source crossed a microcaustic) could be a structure mainly caused by 
a single straight fold caustic, i.e., the curvature and other possible close microcaustics do 
not significantly perturb the simple magnification law. We adopt this last point of view, 
and take the GLITP light curve for image A, including only data points very close to the 
maximum of the HME, to be fitted to the microlensing curves resulting from sources crossing 
a single straight fold caustic. Wc remark that the high asymmetry of the peak as well as 
probabilistic arguments arc two strong reasons against an interpretation of the microlensing 
peak based on a source passing close to a single cusp caustic. 

In Section 2 we present the expected microlensing light curves when axially symmetric 
sources cross a single straight fold caustic, and the fitting procedure. We use a set of 
axisymmetric sources: brightness distributions enhanced at the centre of the source (standard 
physical profile, Gaussian profile, and p — 3/2, 5/2 power-law profiles) and the uniform 
brightness distribution (e.g., Shakura & Sunyaev 1973; Schneider & Weiss 1987; Shalyapin 
2001). Section 3 is devoted to the parameter estimation from the comparison between the 
GLITP microlensing peak in the component A and the expected time evolutions for the 
different source models. A discussion on the V-band and R-band source sizes, the source 
size ratio {Rv/Rr), and the reliability of the source models is also included in Section 3. In 
this paper, to obtain information on the dimension of the V^-band and i?-band sources, we 
will use the measurements of the transverse galactic velocity reported by Wyithe, Webster 
& Turner (1999). Finally, in Section 4 we summarize our results and conclusions. 



2. Theoretical microlensing light curves and fitting procedure 

When a source crosses a single straight fold caustic, a microlensing violent phenomenon 
occurs. This is seen as an important fluctuation in the flux of an image of the source. We 
are going to describe the properties of a family of axisymmetric sources, present the time 
evolution of the fiux (microlensing curve) corresponding to each source model, and finally, 
introduce the methodology to compare the theoretical curves with observational data. 
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2.1. Source models 

2.1.1. Uniform and Gaussian disks 

In a given optical band (C/, B, V, R, I, ...), the two simplest surface brightness dis- 
tributions are a uniform disk and a Gaussian disk (e.g., Schneider & Weiss 1987). For the 
source of uniform intensity, one has a well-defined source radius Ropt{source) that coin- 
cides with the typical radius of the intensity distribution Ropt- Both radii contain the total 
brightness of the source. However, for the source with Gaussian intensity distribution and 
other sources with distributions enhanced at the central regions, we know a typical radius 
of the 2D profile (Ropt), and assume that Ropt{source) ^ Ropt- This last hypothesis per- 
mits us to justify the approximation ~ J^°°_ Pq^. example, in the Gaussian case 
with /opt(r) = e/optexp(— r^/i?^pj, the total brightness included in the circle with radius 
Ropt{source) is given by Iopt{source) — Ioo{l — exp[— i?Qp((soiirce)/i?opJ}, where /oo is the 
brightness derived from an integration between r — and r — oo. If Ropt{source) ^ Ropt, we 
have Iopt{source) ^ I^o- We note that I^^pt = Iopt{Ropt) is the intensity at the typical radius 
Ropt (i.e., a typical intensity), and so, /opt(O) = elgpt for the Gaussian law. As different mod- 
els have typical radii with different meanings, in a unified scheme, we define two source sizes: 
the source radius containing half of the total brightness, Ropt{50%), and the source radius 
containing 90% of the total brightness, Ropt{90%). Using factors /c(50%) and A;(90%), which 
depend on the source model, the radii i?opt(50%) and Ropt{90%) can be expressed in terms 
of the typical radius Ropt [e.g., -Ropt(50%) = k{50%)Ropt\- For a uniform disk, the factors are 
A;(50%) = 1/V2 and A;(90%) = 3/VTO, while for a Gaussian disk, we obtain A;(50%) = Vh[2 
and A; (90%) = VhOO. 

2.1.2. Power-law models 
A more interesting set of source models is the family of intensity profiles with behaviours 

(Popt > 1) 

loptir) = 2^'-*7„,,(l + ryR%r"'^K (1) 

Along with the simplicity of the trends, these power-law models allow to calculate the ex- 
pected microlensing curves in an analytical way (see Shalyapin 2001, and the next subsection 
2.2.). We shall use two particular values of the power index, popt — 3/2, 5/2, which lead to 
specially simple laws for the change in flux during a caustic crossing. 

For these power-law models, the brightness enclosed within the circle with radius r will 
be Iopt{< r) = 7oo[l - (1 + ryRlpty-P°-*], being « Iopt{source). When r = Ropt{50%) 
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the second term in the brackets must be equal to 1/2. So k{50%) = [2^/^p°p* ^) - l]^/^. On 
the other hand, it is easy to show that A;(90%) = [loV(Pcpt-i) _ i]V2. 



2.1.3. Standard accretion disk 

Up to now we described models that are not directly related to physical ideas about the 
quasar central engine, i.e., we dealt with effective models. However, we can also consider 
the standard model of an accretion disk around a supermassive black hole. The standard 
Newtonian model (Shakura & Sunyaev 1973) is based on the supposition that the released 
gravitational energy is emitted as a multitemperature blackbody radiation. The temperature 
profile is 



Ts{r) 



3 GM ■ f jVin^^^'^ 

-Mil - X — 




(2) 



where a is the Stefan constant, G is the gravitation constant, M is the mass of the central 
black hole, and M is the accretion rate. Here, r^n is the inner radius of the accretion disk, 
which is usually assumed to be thrice the Schwarzschild radius of the black hole (rschw)- The 
emitted intensity obeys a Planck law 



and taking into account the redshift of the source Zg, the observed intensity / — Is{l + z, 
is also a Planck function at the temperature T = Tg/ (1 + Zs)i 

2hu^ 1 

^ '~ e^^{hv/kT)-l' ^' 

At a frequency Vopt (the central frequency of the optical filter), from Eqs. (2) and (4) we 
infer 

lovtir) = , (e--l)V 

exp{(r/7?,,,)3/4 [1 - u{r/Ro,t)-'/T ' } " 1 

where V = hptiRopt) = ^hvlpt/i^^"" ' l)c', Ropt = {3GMM /8nay/^[k/huopt{l + Zs)f/^ 
u = {rin/ Ropty^'^, and w = {I — u)~^^*. The parameter u is taken for simplicity to be 
negligible, i.e., -u = (w = 1), in such a way that the final intensity profile leads to a 
theoretical microlensing variation quite similar to the exact standard variation at m << 1. 
At u ~ 1 the simplified version of the standard profile is not a useful approach, but as it will 
be discussed at the end of Section 3, the observed hght curves do not show the strong break 
predicted by an exact model with ~ 1 and the simplified profile {u — 0) works reasonably 
well. 
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For the standard source, the coefficients A;(50%) and A;(90%) must be estimated in a 
numerical way. Their values are: A;(50%) = 2.386, A;(90%) = 7.038. 



2.2. Microlensing curves during a caustic crossing 

The microlensing light curves are basically calculated by convolving intensity distribu- 
tions with the magnification pattern associated with a single straight fold caustic. This 
simple magnification pattern is given by (e.g., Schneider & Weiss 1987): a constant back- 
ground magnification Aq at points located outside the caustic, and a law Aq + ac/Vd at 
points that are placed inside the caustic. Here, qq is the caustic strength and d is the 
perpendicular distance to the caustic line. 

We take a coordinate frame in which the caustic line is defined by the y-axis and 

the centre of the circular source has the coordinates {xc,0). To obtain the flux from a 
part of the disk at {x,y), its intrinsic intensity should be multiplied by both the extinction 
factor eopt and the solid angle dQ it subtends on the sky. Considering the relationship: 
dQ = A{x,y)dQ^, where A{x,y) = Aq + acH{x)/ ^/x is the achromatic magnification factor, 
H{x) is the Heavisidc step function, and dQ^: is the solid angle in the absence of lens, the 
elemental fiux will be eoptA{x,y)Iopt{r)d^}^,, with r = \/ {x — Xc)"^ + y"^- Integrating over 
the sky, it is inferred a radiation fiux of the QSO image (e.g., Jaroszyhski, Wambsganss & 
Paczyhski 1992) 

Fopt{xc) = ^ j j A{x,y)Iopt{r)dxdy, (6) 

where dQ^ = dxdy/D^ and Dg is the angular diameter distance to the source. This fiux can 
be rewritten as 

ac 



^optyXc) — 



1 + —r-fopt{xc) 

^0 



(7) 



where Qopt — J J Iopt{''^)dxdy represents the total intrinsic brightness of the source and 

, . _ J J H{x)x-^/^Ioptir)dxdy 
J J Ioptir)dxdy ■ 

Finally, using normalized coordinates ^ = x/Ropt and rj = y/Ropt, one finds gopt = Rlpt^ptK 
and fopt{xc) = R^/'^J{xc/Ropt), being 

K^^J j loptiRopt^ (e - xjRoptf + v'Hdr), (9) 

and 

J J Iopt[RoptV (C - Xc/Roptf + V^]d^dr] 
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The function J{xc/Ropt) was also used by Schneider & Weiss (1987), Shalyapin (2001), and 
other authors. 

Prom the results in the previous paragraph, one obtains 

F^t{xc)^Fo + FcJ{z), (11) 

where Fq = {Rl^J DDeoptAoIoptK , Fc = (Fo/y^)(acMo), and z = Xc/Ropt- In Eq. 
(11) two chromatic amplitudes appear: the background flux (Fq) and the amplitude of the 
contribution due to the extra magnification inside the caustic {Fc). On the other hand, the 

constant K and the function J{z) depend on the source model, and we focused on the shape 
factor J . In particular we estimated the behaviours of J{z) corresponding to the five models 
discussed in the previous subsection. For the uniform, Gaussian and standard disks, J{z) 
was deduced in a numerical way, whereas for the Popt = 3/2, 5/2 power-law profiles, J{z) 
has an analytical form. In the case popt = 3/2, Shalyapin (2001) gave a simple expression 
for the shape factor, and when p^pt = 5/2, 



Jiz) = ^^^ + ZZ?- , . (12) 

^ ' [2(l + ^2)(^l^_^)]3/2 V ) 

The trends of J{z) for the five models are depicted in Fig. 1: uniform disk (dotted fine), 
Gaussian disk (dashed line), Popt — 5/2 power- law model (dash-dotted line), Popt — 3/2 
power-law model (dash-three-dotted line), and standard accretion disk (solid line). To derive 
microlensing light curves, the final step will be to use the trajectory of the centre of the source: 
Xc{t) = V±{t — to), where V± is the quasar velocity perpendicular to the caustic line, and to 
is the time of caustic crossing by the source centre. We implicitly assumed that the source 
enters the caustic, i.e., > at t > to- Moreover, we remark that the time t is measured 
by the observer. Inserting the trajectory into Eq. (11), we infer microlensing curves 

Fopt{t)^Fo + Fcj(^-^y (13) 

where At = Ropt/V± is the crossing time. When the main part of an axisymmetric source 
(i.e., the circle containing half of the total brightness) crosses a fold caustic, Eq. (13) will 
give a good approach to the flux of the involved image as a function of time. However, 
the photometric behaviour relatively far from a particular fold caustic could be different to 
the law (13). Firstly, the simple magnification factor A{x,y) = Aq + acH{x)/ y/x may be 
perturbed as due to the curvature of the fold caustic. This problem was studied by Fluke 
& Webster (1999), and more recently, by Gaudi & Fetters (2002). Secondly, the presence 
of another caustic (i.e., the existence of a network of caustics) could dramaticly change the 
magnification pattern A{x, y) and the predicted time evolution of the flux. Thirdly, the 
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background magnification ^40, which is caused by all the point source images not associated 
with the fold of interest, may actually be a function of the position: Aq = Aq{x, y). Gaudi 
& Fetters (2002) introduced microlensing curves incorporating a slowly varying background. 
In practice, when one deals with a given observed microlensing peak or event, the trends (13) 
must be the first choices to fit it. Of course a family of reasonable source models should be 
tested. After to do the initial fits, if there is evidence for important post-fit residues using all 
the effective and physical scenarios, then some corrections must be made (curvature, another 
caustic, and so on). As we will see in Section 3, our dataset does not suggest the need of 
corrections, and the behaviour (13) works very well when some source models are considered. 



2.3. Theory vs. observations 

Any microlensing peak or event in a component of a lensed quasar can be compared 
with the theoretical hght curves presented here above [see Eq. (13)]. Given a source model, 
the light curve during a caustic crossing (isolated and straight fold caustic) depends on 4 
parameters: 

1. Fo - a background flux (e.g., in mJy); 

2. Fc - Si flux related to the extra magniflcation inside the caustic (e.g., in mJy); 

3. to - a, time of caustic crossing by the source centre (e.g., in JD-2450000); 

4. At - a. typical crossing time, which is deflned as the ratio between the typical radius 
Ropt and the perpendicular motion V± (e.g., in days). 

So, assuming a particular intensity proflle, the task is to estimate the values of the pa- 
rameters Fq, Fc, to, and At which best describe the observed behaviour. The estimation 
of model parameters is carried out from a fitting method. To fit N observational data 
-fopt(l)i •••) Fopt{N) with errors ai, a^, respectively, to the expected ones at times ti, .... t]v- 
-Fo + FcJ[{ti — to)/ At], 1 < i < N , a, chi-square minimization will be used. We arc going to 
search for the values of the four free parameters {Fo, Fc, to. At) which minimize the sum 

X-(F„. Fa. t,, M)^t{ ^^'t') - ^° - - r . (14) 

i=l ^ J 



We must fit a theoretical law that is linear in two parameters {Fo, Fc) and non- linear in 
to and At. On the other hand, to flnd the best values of Fo, Fc, to, and At, one must solve 
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the system of equations dy^ jdF^ — dy^/dFc — jdi^ — dy^ /dlH^t — 0. As the function 
Fopt{t) is hnear in Fq and Fc, the equations dx^/dFQ = dx^/dFc = lead to analytic 
relations: Fq = Fo{to, At) and Fc — Fc{to,A.t). From Eq. (14) and these constraints, we 
can make a new chi-square 



It is a clear matter that the function xt is the old function which has been minimized 
with respect to the two linear parameters, and thus, the complex problem involving 4 free 
parameters is reduced to a very simple problem: a minimization of the 2-dimensional dis- 
tribution x^(to,At). Assuming we have the best values for {Fq, Fc,to, At), it is necessary 
to estimate the errors on the parameters. In particular, we concentrated on the most rele- 
vant parameter At. Given the chi-square global minimum x^(min) = x2(min), to infer the 
confidence intervals of a single parameter At, we take the different values of the parameter 
of interest verifying the conditions Ax^ — xtitojAt) — x^(min) < 1 (68 percent confidence 
interval, i.e., la), Ax^ < 4 (95 percent confidence interval, i.e., 2a), and so on. 

To complete the process, we must have an idea of the quality of the fit. If the data 
correspond to the theoretical law and the deviations (due to the observational noise) are 
Gaussian, x^ should be expected to follow a chi-square distribution with mean value equal 
to the degrees of freedom, dof = N — 4. We thus expect to be closed to iV — 4 if the fit is 
good. A quick test is to form the reduced chi-square x^ = x^/dof, which must be close to 1 
for a good fit. Moreover, for dof > 30, the chi-square distribution is essentially normal with 
standard deviation of \/2 dof. Therefore, the relative deviation S = \x^ — dof\j\/2 dof is 
expected to be < 1. 

3. Confrontation between GLITP data for Q2237-K0305A and theoretical light 



The results of the comparison between the GLITP y-band light curve for Q2237+0305A 
and the five expected time evolutions are presented in Table 1, while the results from the 
comparison in the R band appear in Table 2. We are mainly interested in measurements of the 
y-band and it!-band source sizes, and so, only the uncertainties {la intervals) in the estimates 
of the relevant parameter At are quoted. In Table 3 (second column) we have included the la 
confidence intervals of the source size ratio (g = Rv/Rr). Table 3 also includes upper limits 
on the half-light radii. As Ropt{50%) = A;(50%)Vj_At, to infer the constraints, we used V± < 
7000 km s"^ (Wyithe, Webster & Turner 1999) and At < Atmax = Attest- fit + (7+(Ai). For 




(15) 



curves 
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some source models, information on the source radii (containing 100% of the total brightness) 
or Ropt(90%) in the V and R bands, is quoted for discussion. Wyithe, Webster & Turner 
(1999) determined upper limits on the transverse galactic velocity {vt) from the distribution 
of light-curve derivatives. In Fig. 12 (bottom panels) of Wyithe, Webster & Turner, it was 
presented the 95 per cent upper limit to Vt as a function of the adopted random errors, 
the microlens mass distribution, and the source details. For mass distributions which are 
characterized by a mean microlens mass equal to or less than 1 Mq, and any reasonable 
choice of the random uncertainties, the global bound is of Vt < 900 km s~^. This bound is 
valid for any source size and intensity profile, and we assumed it as an absolute upper limit. 
Using the relationship V — {Ds/D^jVt, where is the angular diameter distance to the 
lens (deflector) and V is the transverse quasar velocity, one easily derives that Vj^ < V < 
7000 km s-i (0=1). 

The upper hmits on the F-band and it!-band source sizes are presented in Table 3. For 
the uniform disk: Rv,Rr < 6.3 10"^ pc (= 1.9 10^^ cm) . For the p = 3/2 power- law model 
and the standard model, the radii containing 50% of the total disk brightness may be as 
large as 2.7 10^^ pc, and interestingly enough the radii containing almost all the brightness 
(90%) must be smaller than 1.2 10~^ pc (= 3.7 10^^ cm), in a reasonable agreement with the 
expected radial size of a hot accretion disk around a 10*^ Mq black hole (10^ rschw ~ 10~^ 
pc) . With respect to the parameter q, a value of q less than 1 would suggest the existence of 
two concentric sources with different radii, the y-band source being inside the it!-band one. 
Unfortunately, although the preferred values of q vary in the interval 0.65 - 0.98, i.e., q < 1, 
the error bars are quite large and we cannot fairly distinguish between the cases q < 1 and q 
= 1. Only the Gaussian model led to g < 0.86 at la confidence level. For the standard disk, 
Ropt oc i^opt^ Q = (^i?/^v^)^^^ ~ 0.8. However, this standard source size ratio cannot be 
confirmed from our indirect measurement. 

In Table 1, the 5{min) values are not good for the y-band uniform and Gaussian disks 
(see subsection 2.3). However, it is evident that both the pv = 3/2 power-law model and 
the V-band standard disk work very well. The intensity distribution with a power-index of 
3/2 and the standard accretion disk are also favored from data in the R band (see Table 
2). For the standard source, curves of constant chi-square in the to (JD-2450000)-At (days) 
plane are showed in Fig. 2. In Fig. 2 we can see the contours associated with Ax^ = 
1 and Ax^ = 4. On the other hand, in Fig. 3 we drawn together the observed hght 
curves and the corresponding standard fits. The agreement between GLITP observations 
and fits is excellent, and taking as reference the day 1500, it is evident the existence of 
an important asymmetry. For comparison, in the top panel (dashed line), we also drawn 
the best-fit from the theoretical microlensing curve related to the exact standard accretion 
disk (including realistic edges). As the exact microlensing law fits the observations slightly 
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better than the approximated one, we confirm the high feasibihty of the model based on 
physical grounds. However, although the standard disk led to values of x^lmin) ~ 1, the p 
= 3/2 power-law model is also in good agreement with the observations. As a consequence 
of this result, we tested if the p = 3/2 power-law model is a simplified version of the exact 
standard model or, on the contrary, we are handling two very different scenarios. The most 
simple test is a comparison of the ID intensity profiles, i.e., the integrals of the 2D intensity 
distributions along the direction defined by the caustic hne. Therefore we first simulated a 
2D intensity distribution associated with a standard accretion disk (exact model) . Then they 
were derived the standard ID profile and all the effective ID profiles fitting the standard 
one. The comparison between the standard behaviour and the effective trends appears in 
Fig. 4. The filled circles represent the standard ID profile (including the two typical "horns" 
close to a; = 0), while the solid lines trace the power-law {p = 3/2) ID profile (left-hand top 
panel), the power-law (p = 5/2) ID profile (right-hand top panel), the Gaussian ID profile 
(left-hand bottom panel), and the uniform ID profile (right-hand bottom panel). It is now 
clear that the p = 3/2 power- law model very closely mimics the behaviour of the standard 
scenario, i.e., it is a rough version of the exact standard model. 



4. Conclusions 

We have analyzed the peak of a microlensing high-magnification event which was accu- 
rately monitored by the GLITP collaboration (Alcalde et al. 2002). The prominent event 
has occurred in image A of the quadruple system QSO 2237-1-0305 (Wozniak et al. 2000b), 
and the GLITP team observed its peak in two optical bands {V and R). As both F-band 
and it!-band light curves are characterized by high flux and time resolutions, we attempted 
to interpret the observational trends in terms of some theoretical microlensing light curves. 
In principle, the observed microlensing peak could mainly result from the source crossing a 
microcaustic or passing close to a microcusp. However, taking into account probabilistic ar- 
guments and the important asymmetry around day 1500 (see Fig. 3), we have discarded the 
hypothesis of a source in the vicinity of a single cusp caustic. Thus we only considered the 
most simple alternative picture: a source crossing a single straight fold caustic. Different ax- 
ially symmetric source models were used, and consequently, several theoretical microlensing 
curves were compared with the observed ones. To discuss the reliability/feasibility of differ- 
ent intrinsic intensity proflles, all source models were chosen to cause theoretical microlensing 
curves with the same number of free parameters. These parameters are two characteristic 
fluxes, the time of caustic crossing by the source centre, and the typical crossing time. 

Our main results and conclusions are: 
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1. Prom the fits and the upper hmits on the transverse galactic velocity claimed by 
Wyithe, Webster & Turner (1999), we inferred constraints on the K-band and i?-band 
source sizes. For a uniform source model, the V"-band and -R-band radii should be less 
than 6.3 10~^ pc = 1.9 10^^ cm, at la confidence level. For a Gaussian disk, the typical 
y-band radius is < 7.8 10^"' pc (= 2.4 10^^ cm, at la CL), while the typical i?-band 
radius can be a little larger (< 1.6 10^^ pc = 4.9 10^^ cm, at la CL). Taking the total 
size of the i?-band source as the source diameter for a uniform disk, or the full-width 
at one-tenth maximum for a Gaussian profile, we found that the la upper limits on 
the i?-band source size are 3.9 10^^ cm (top-hat) and 1.5 10^^ cm (Gaussian). These 
bounds agree well with the results by Wyithe et al. (2000). On the other hand, for the 
standard accretion disk, we obtained that 90% of the K-band and i?-band luminosities 
are radiated within a radius of 8 10^"^ pc (= 2.5 10^^ cm, at la CL), 1.2 10~^ pc (= 
3.7 10^^ cm, at 2a CL). As the total luminosity of a standard accretion disk around 
a 10^ Mq black hole will be probably enclosed in a circle with radial size of ^ 10~^ 
pc, the standard results are highly consistent with the current paradigm on the central 
engine in QSOs. However, we remark that other central masses also agree with the 
constraints. Once upper limits of ~ 10^^ - 10^^ cm are known, we may attempt to test 
the hypothesis of a neghgible caustic curvature (at Vji < 7000 km s~^, the parallel path 
length during ~ 100 days will be also less than 10^'^ - 10^^ cm). The typical caustic 
curvature radius is usually assumed to be the Einstein-ring radius on the source plane, 
i.e., Rcc ~ Re ~ lO^^(m/M0)^/^ cm, where m is the microlens mass = 1, H = 
60 - 70 km s^^ Mpc^^). Therefore, for m ~ IMq, the ratios between the optical 
radii and the typical caustic curvature radius will be smaller than 0.01 - 0.1, and this 
result supports the straight fold caustic approximation. For a small microlens mass 
(m ~ O.IMq), we however have a smaller Re and cannot confirm the weakness of the 
curvature effects. 

2. Wc also studied the source size ratio, i.e., the ratio between the y-band radius and the 
i?-band radius. The source size ratio values are independent of the transverse quasar 
velocity. At la CL, only the Gaussian source model led to a F-band source being 
inside the i?-band one. For the standard source, one hopes for a ratio of about 0.8, but 
this expected value could not be confirmed from our microlensing experiment. The la 
confidence interval (0.53 - 1.26) is in agreement with the two possible situations: the 
\^-band source being inside the R-hand source, and both sources having a similar size. 

3. An important issue is the reliability /feasibility of several source models tested by us. 
With respect to this point, we deduced very interesting conclusions. The usual top-hat 
and Gaussian profiles are not favored from the data in the V and R bands. The results 
are better when power-law profiles are assumed, particularly for the model with a power 
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index of 1.5, which more closely resembles the one-dimensional intensity profile of the 
exact standard accretion disk. For the standard source model, the reduced chi-square 
values are very close to one. So, an accretion disk around a supermassive black hole 
seems a good candidate to be the optical continuum main source in QSO 2237+0305 
(see Fig. 3), and a measurement of the black hole mass and the mass accretion rate 
could be made in that system (e.g., Yonehara et al. 1998). This last topic will be 
discussed in a separate paper. We also note that a hybrid scenario in which the light 
in the V band is emitted from an accretion disk and the i?-band light comes from the 
accretion disk and another extended region, is not in contradiction with the observed 
microlensing peaks. This result totally agrees with the main conclusion of Jaroszyhski, 
Wambsganss & Paczyhski (1992), who previously tested the hybrid model from old 
observational data. The assumption of a complementary i?-band extended source only 
modifies the expression of the i?-band background fiux Fqr in Eqs. (11) and (13), and 
thus, the interpretation of its best-value and uncertainties. Exclusively using GLITP 
data, one can try to analyze the possible existence of a light excess in the R band. To 
do this task, it is needful an accurate calibration of the flux in both optical bands and 
small uncertainties in the measurements of the background fluxes. A detailed study 
of this topic is however out of the scope of the paper, which is focused on the optical 
continuum main (compact) source. On the other hand, in order to flt the observed 
spectrum of the quasar. Ranch & Blandford (1991) did not consider additional R 
light coming from a large region, but adopted a non-classical accretion disk model. 
Unfortunately, the disk model by Ranch & Blandford (1991) cannot account for the 
old microlensing variations in QSO 2237-1-0305. 

We remark that the OGLE collaboration also monitored the F-band event in Q2237-I-0305A 
(http://www.astro.princeton.edu/~ogle/ogle2/huchra.html). The GLITP F-band photom- 
etry traced the peak of the microlensing event, whereas the OGLE F-band dataset described 
the behaviour of the whole fluctuation. In comparison with the OGLE observational proce- 
dure, the GLITP observations are of higher quality because they were obtained with a larger 
telescope, using a detector with better resolution, and on nights with better seeing. There- 
fore, as due to the quality of the observations and the excellent sampling rate, the GLITP 
V-band peak is probably the best tracer of the underlying signal around the maximum of the 
whole flux variation. On the other hand, the whole event seems to be caused by a complex 
magniflcation, which may include ingredients such as the curvature of the fold caustic (e.g.. 
Fluke & Webster 1999), the presence of another caustic, and a non-constant background 
magniflcation (e.g., Gaudi & Peters 2002). The simple flt to the whole event may thus give 
a wrong estimation of the parameters. In any case, we chose two OGLE observation periods 
to infer standard solutions and compare them with the standard flt presented in Table 1. 
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The first period included only the end of 1999, from day 1450.6 to day 1529.5. The dataset 
is called OGLE99, and it corresponds to the GLITP monitoring period. The second dataset 
(OGLE99-00) covers the 1999-2000 seasons, more exactly from day 1289.9 to day 1766.7. 
In Table 4 we can see all at once the fits from the GLITP, OGLE99, and OGLE99-00 light 
curves. As x^{min) > 1 for the OGLE best solutions, in the time parameter estimation from 
the OGLE datasets, we considered the la bounds associated with Ax^ = x^{min) rather 
than Ax^ = 1 (e.g., Grogin & Narayan 1996). The GLITP and OGLE99 fits are not consis- 
tent each other, but the amplitudes and the crossing time from the OGLE99-00 brightness 
record are close to the values of Fq, Fc, and At from the GLITP dataset. 

Finally, we must emphasize that an important progress can be made from accurate and 
detailed data of a microlensing peak. One can discuss on the reliability/ feasibility of different 
source models leading to the same number of free parameters, and thus, to discard some of 
them and to find the models that agree with the observations. Given a good model, which 
is consistent with the data, it is possible to obtain a robust upper limit on the size of the 
optical source. If, for example, the good model is an accretion disk around a massive black 
hole, then one can try to work out a technique to measure the central mass and the accretion 
rate. Moreover, as the accurate and well-sampled microlensing peaks seem to be inconsistent 
with some intensity profiles, it could be very interesting to apply the deconvolution method 
to these peaks and to derive the best brightness distribution in a more direct way (Grieger, 
Kayser & Schramm 1991; Agol & Krohk 1999; Mineshige & Yonehara 1999). Although we 
had not success in the accurate and robust indirect estimation of the ratio between the V- 
band radius and the it!-band radius (the parameter q) , new multiband monitoring campaigns 
could also lead to relevant measurements of source size ratios. With regard to this last issue, 
we note that a direct estimate of q from the cross-correlation of our \/-band and i?-band 
light curves is now in progress, but the expectations are not very promising. 
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FIGURE CAPTIONS 

Fig. 1. — Shape factor J{z). Five source models are considered: uniform disk (dotted line), 
Gaussian disk (dashed line), Popt = 5/2 power- law model (dash-dotted line), Popt = 3/2 
power-law model (dash-three-dotted line), and standard accretion disk (solid line). 

Fig. 2. — Ax^ = 1 (interior solid lines) and Ax^ = 4 (exterior solid lines) contours in the 
to (JD-2450000)-At (days) plane. The source is assumed to be a standard accretion disk. 
The datasets are: the GLITP V-hand light curve of Q2237+0305A (top), and the GLITP 
i?-band light curve of Q2237+0305A (bottom). 

Fig. 3. — Observed light curves and the corresponding standard fits. Top: y-band. Bottom: 
i?-band. The dashed line in the top panel (V^-band) represents the best-fit from an exact 
standard accretion disk crossing a caustic line, and the weak break around day 1480 is an 
edge effect. We can see a good agreement between both V-band fits. 

Fig. 4. — Comparison of a given exact standard ID intensity profile (filled circles) and the 
effective ID profiles fitting it (solid lines). The effective source models are: p = 3/2 power- 
law (left-hand top panel), p = 5/2 power-law (right-hand top panel), Gaussian (left-hand 
bottom panel), and top-hat (right-hand bottom panel). 
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Table 1: Fits to the GLITP V-hand light curve for Q2237+0305A. 



Source model 


Fo (mJy) 


Fc (mJy) 


to (JD-2450000) 


At (days) 


X^{min)/6{min) 


Uniform 




0.72 


0.07 


1489.6 


30.2; 


fl.3 
-1.3 


1.33/1.62 


Gaussian 




0.71 


0.09 


1488.6 


31.0; 


f8.3 
-2.8 


1.31/1.52 


Power-law {pv = 


5/2) 


0.68 


0.11 


1486.1 




1.18/0.88 


Power-law {pv — 


3/2) 


0.65 


0.19 


1484.7 


33.1; 


f8.3 
-6.1 


1.08/0.39 


Standard 




0.59 


0.30 


1480.6 


39.6i^°5" 


0.99/0.05 
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Table 2: Fits to the GLITP i?-band light curve for Q2237+0305Af 



Source model Fq (mJy) Fc (mJy) to (JD-2450000) At (days) (min) / 6 {min) 



Uniform 




1.28 


0.09 


1489.3 




1.44/2.09 


Gaussian 




1.18 


0.17 


1483.7 


47.61E1,^ 


1.30/1.42 


Power-law {pr — 


:5/2) 


1.14 


0.20 


1481.6 


a,r. r.+UA 
00.0_j8 


1.19/0.90 


Power-law {pr — 


:3/2) 


1.12 


0.33 


1481.4 




1.12/0.57 


Standard 




1.06 


0.46 


1478.7 




0.99/0.05 



°Thc observed fliix at day 1473 deviates from neighbouring points and it was removed from our brightness 
record 



-27- 



Tablc 3: Source size ratios and boLUids on the (limcusiou uf'tlic 1 -baud and /iL-baii(l sources. 



Source model 


q — Rv/ Rr 


i?y(50%) X 10^ (pc) 


i?i?(50%) X 10^ (pc) 


Uniform 




n 98+°-°^ 


< 0.044 (0.061)^ 


< 0.044 (0.063)"^ 


Gaussian 




U.UO_Q 47 


< 0.065 


< 0.133 


Power-law {p = 


5/2) 




< 0.091 


< 0.150 


Power-law {p = 


3/2) 




< 0.140 (0.807)^ 


< 0.205 (1.172)^^ 


Standard 






< 0.234 (0.689)^ 


< 0.271 (0.800)^^ 



"Upper limit on Ry x 10^ (pc) 
''Upper limit on Rr x 10^ (pc) 
'^Upper limit on i?y(90%) x 10^ (pc) 
''Upper limit on i?fl(90%) x 10^ (pc) 
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Table 4: Comparison between the standard fits to the GLITP and OGLE data. 



Dataset 


Fo (mJy) 


Fc (mJy) 


to (JD-2450000) 


At (days) 


X^{min) 


GLITP 


0.59 


0.30 




39.6+^'^/ 


0.99 


OGLE99 


0.72 


0.20 


1487.61}:^ 


14.311^ 


L51 


OGLE99-00 


0.60 


0.34 


1471.211^ 




7.49 



